Charge degree of freedom and single-spin fluid model in YBa2Cu40s 
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We present a 17 nuclear magnetic resonance study in the stoichiometric superconductor 
YBa2Cu408- A double irradiation method enables us to show that, below around 200 K, the 
spin-lattice relaxation rate of plane oxygen is not only driven by magnetic, but also significantly 
by quadrupolar fluctuations, i.e. low-frequency charge fluctuations. In the superconducting state, 
on lowering the temperature, the quadrupolar relaxation diminishes faster than the magnetic one. 
These findings show that, with the opening of the pseudo spin gap, a charge degree of freedom of 
mainly oxygen character is present in the electronic low-energy excitation spectrum. 



PACS numbers: 7I.45.-d, 74.25.Nf, 74.72.Bk, 76.60, 

One of the central issues in understanding cuprate- 
based high-temperature superconductors (HTSC) has 
been to determine the minimal number of electronic de- 
grees of freedom which are necessary to describe the 
physics at the atomic scale of these structures. Zhang and 
Rice proposed that doped holes, which go into oxygen 
states, form a spin-resonant singlet state with the quasi- 
localized holes at copper Therefore, only one spin 
degree of freedom would be necessary to describe the low- 
energy dynamics of the electronic system in the normal 
state of the HTSC. This scenario is often called single- 
spin fluid (SSF) and is described within a t-J model. 
However, as was discussed by various authors [0, this 
model might be an oversimplification. 

Due to its local character, nuclear magnetic resonance 
(NMR) is a very attractive method to address these 
questions. NMR probes the low-energy excitations of 
the HTSC electronic system and it provides information 
about copper and oxygen independently. Studies of the 
uniform spin susceptibility, Re{x(<Z = 0, w — 0} |T|], indi- 
cate that only one spin degree of freedom is involved thus 
supporting the SSF model. On the other hand, oxygen 
and copper reveal a very different temperature depen- 
dence of the spin-lattice relaxation B which is related to 
the dynamic spin susceptibility, Im{x(<7, Wl)} ||. This 
could suggest independent spin degrees of freedom. How- 
ever, as pointed out before the hyperfine field, due to 
antifcrromagnetically correlated copper spins, cancels at 
the oxygen site and this fact can explain the different 
temperature dependences within a SSF model. 

There is, however, evidence that this view is incom- 
plete. The 17 NMR in Y-Ba-Cu-0 HTSC has exten- 
sively been studied by various groups. Especially the 
pronounced temperature dependence of the anisotropy 
of the experimentally determined effective spin-lattice re- 
laxation rate, 17 W c g, of the plane oxygen is still puzzling 
while the SSF model predicts an almost tempera- 
ture independent rate anisotropy. (The meaning of effec- 
tive will be discussed below.) As we pointed out in Ref. 
p|, it is also very difficult to explain, within a simple 
SSF model, the ratio of this rate and the yttrium rate, 
namely 17 W e s / 89 W (with the external magnetic field Bq 



parallel to the c axis). Furthermore, in a detailed anal- 
ysis, Walstedt et al. [[Io| compared their 17 and 63 Cu 
NMR results from La2- x Sr x Cu04 studies with inelastic 
neutron scattering data and concluded that a single-band 
picture is inadequate. 

In this Letter, we will show that the plane oxygen spin- 
lattice relaxation in Yb^CuziOs, below about 200 K, 
is not only driven by magnetic but also by quadrupolar 
fluctuations, i.e. low-frequency charge fluctuations. This 
fact implies that in the temperature region where the 
pseudo spin gap is present, a charge degree of freedom has 
to be taken into account, in addition to the spin degree 
of freedom; hence the simple SSF model is inadequate. 

Our results were obtained by a double irradiation 
technique which we developed to extract the quadrupo- 
lar contribution to the overall NMR relaxation O]; we 
sketch this technique as follows. The relaxation of the 
spin system towards its thermodynamic equilibrium is 
described by the so-called master equation |L2j and, in 
the interaction representation |R||l4| , it is given by 



R(W,Wi,W2)[P(t)-P(0)] + SrfP(t). (f) 



Here, P(t) is the nuclear spin system population vector 
of the different energy levels with P(0) being the equi- 
librium value. The relaxation matrix R(Wi W\,Wz) is, 
in general, a function of the magnetic relaxation rate 
W (causing transitions with |Am| = I) and the two 
quadrupolar relaxation rates W\ (|Am| = I) and W% 
(| Am | = 2) |T^,|I^]. The transition probabilities for a 
spin I = 5/2 (like 17 0) are sketched in Fig. S rf de- 
notes transitions due to an additional stimulating radio- 
frequency (rf) field which is essential for the separation 
of magnetic from quadrupolar contributions to the relax- 
ation, as will be discussed below. In HTSC, the Zeeman 
degeneracy is lifted due to the crystal field so that differ- 
ent transitions have distinguishable resonance frequencies 
[jl2|| . In case the relaxation is purely magnetic it is easy 
to show that the saturation of any line does not affect the 
intensities of others. This is no more true for the general 
case where W as well as W\ and W 2 contribute to the 
relaxation. 
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FIG. 1. Probabilities for transitions between the spin 
energy levels effected by magnetic (W) and quadrupolar 
(Wi, W2) spin-lattice relaxation processes for / = 5/2. SW 
denotes transitions due to a stimulating radio- frequency field. 

In the limit of dynamical saturation (i.e. dP(t)/dt — 
0), Eq. (|l|) can be solved exactly The essential re- 
sult is that, mainly due to the Am| = 2 transitions, the 
intensities of different lines change compared to a stan- 
dard detection method (e.g. spin-echo sequence) where 
5 r f is zero. In our present experiment, Srf saturates the 
(—1/2,1/2) transition and we observe the other transi- 
tions, yielding signals with intensities /dyn- The spin- 
echo sequence (with S r f=0) provides signals with inten- 
sities I se . For the (—3/2,-1/2) transition the enhance- 
ment factor, E, is then given by [O: 



urated the outer high-frequency satellite ]TT[ |. 
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with \x = A5W 2 (WW + 2W X + W 2 ) and C = 400fW 2 
+ lOOOVFPFi + AQWl + H00WW 2 + I6OWW2 + 45W|. 
E is a rather insensitive function with respect to W±, 
since W\ connects the same energy levels as W, except 
the (—1/2, 1/2) transition [|l~Tf| . However, E depends sig- 
nificantly on W 2 which makes the separate detection of 
quadrupole relaxation possible. The enhancement E for 
various transitions is quite different; for the (—5/2, —3/2) 
transition, there is almost no enhancement. More infor- 
mation and details concerning the method and its exper- 
imental realization as well as different cross-checks of the 
results are given in Ref. [jll| . 

The measurements were performed on an oriented and 
17 enriched YBa2Cu40s powder sample, used in previ- 
ous studies in a field of 8.9945 T applied along the c 
axis; Fig. || presents the results for T = 95 K. The chain 
sites 0(1) are not of interest in this work and therefore 
will not be discussed. Also, we are not concerned with 
the plane oxygen satellites 0(2,3) splitting which is due 
to the orthorhombic symmetry |l5| of the crystal. Strik- 
ing the eye is the negative intensity of the central transi- 
tion ( — 1/2, 1/2) in the difference spectrum, which means 
saturation of this transition. This is in contrast with the 
central and the satellite transitions of chain oxygen which 
are not saturated because 0(1) relax considerably faster 
than plane nuclei for which the pulse sequence was opti- 
mized. That this explanation is correct has been proven 
by the symmetric experiment where we dynamically sat- 
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FIG. 2. (a) ir O central transition and high frequency satel- 
lites of chain and plane oxygen in YBa2Cu4 0s, respectively, 
obtained by a standard spin-echo experiment, (b) 17 spec- 
trum obtained after dynamic saturation of the ( — 1/2,1/2) 
transition. Bottom: Difference of spectrum (b) and (a). 

The most important feature of the difference spec- 
trum is the remaining positive intensity at the position of 
the (—3/2, —1/2) transition. This intensity enhancement 
clearly shows that there is a quadrupolar contribution to 
the relaxation of plane oxygen nuclei. At the position of 
the (—5/2, —3/2) transition, the intensity is almost zero, 
as expected from our calculations. 

Fig. H shows, among others, the temperature de- 
pendence of the intensity enhancement, E, of the 
(—3/2,-1/2) transition due to the saturation of the 
(—1/2,1/2) transition. With these values, Eq.(2) al- 
lows us to calculate 17 (W 2 /W). In doing so, we have set 
W\ = W 2 , which is not only reasonable, if one assumes 
more or less isotropic quadrupolar fluctuations, but also 
brings our result in accord with 17 W e g/ S9 W data as we 
will see below. The 17 (W 2 /W) results are plotted in Fig. 
H by using an ordinate axis in such a way that the points 
coincide with the bullets. In oder to estimate 17 W 2 it- 
self, we have re-measured the temperature dependence 
of the nuclear magnetization recovery (obtained by stan- 
dard NMR) and then fitted the data by the theoretical 
expression for combined magnetic and quadrupolar re- 
laxation [|4| using W\ = W 2 and W as parameters and 
our values W 2 /W = W\/W . The results of this fit are 
presented in the inset of Fig. ||. 

We like to stress that 17 W is different from the value 
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of 1 W e g one gets if the fit procedure uses the theoretical 
expression for pure magnetic recovery (Wi = W2 = 0). 
Lacking experimental hints, this fit procedure has been 
applied in all previous HTSC oxygen relaxation studies, 
i.e. quadrupolar oxygen relaxation, if any, has been con- 
sidered to be negligible. Now, having 17 W that prop- 
erly represents the magnetic part of the relaxation, wc 
can evaluate the ratio 17 W/ 89 W and compare it in Fig. 
I with the ratio 17 W cB / 89 W from Ref. §]. Obviously, 
the 17 W/ 89 W temperature dependence becomes much 
weaker than the 17 W cS / 89 W one; it actually approaches 
the theoretical prediction [solid curve in Fig. |] of the 
MMP model Jig) . This agreement is satisfying and it is 
based on the assumption 17 W\ ~ 17 Wi we made above. 
The 17 W2 results have relatively large errors as a conse- 
quence of the large W2/W errors. While the data points 
suggest, above 100 K, an increase of 17 W2 with falling 
temperature, also a temperature independent 17 W2 is 
compatible with the errors. Be that as it may, the es- 
sential result is the presence of 17 W2 rather than its tem- 
perature dependence in the normal state. 
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FIG. 3. Temperature dependence of intensity enhance- 
ment of the 17 (-3/2,-1/2) transition (•); 17 W cft / 89 W 
(«); 17 Wf 9 W (x); theoretical prediction for 17 W / 89 W (solid 
curve) jiff. Inset: 17 W and 17 W2 for plane oxygen. 

What is the origin of the quadrupolar relaxation? 
A simple phononic mechanism can be ruled out for 
two reasons. First, the temperature dependence of the 
quadrupolar contribution is incompatible with a power 
law T a with a = 2 ... 7 |!|L7| expected for relax- 
ation due to phonons. Second, a dramatic softening 
of a phonon mode, that could in principle lead to the 
observed temperature dependence, is not taking place. 
Though, there is a weak softening of phononic modes 
in YBa2Cu40g, as observed by Raman scattering ex- 
periments the effect, besides occuring at a too low 
temperature, is much too small to account for the ob- 



served strong increase of the quadrupolar contribution in 
the oxygen relaxation. We can also exclude quadrupo- 
lar relaxation due to any defect motion. Such a motion 
would be clearly observed in the 137 Ba relaxation since 
the 137 Ba quadrupole moment is much larger than the 
17 one. Yet, we did not detect such an effect [ jl9| . 

A clue to the origin of the quadrupolar relaxation is 
provided by the fact that below about 200 K, if the tem- 
perature decreases towards T c , 17 (W2/W) significantly 
grows and then diminishes in the superconducting state. 
The increase of 17 (W2/W) becomes pronounced around 
a temperature where we recently detected || an 
electronic crossover associated with anomalies in vari- 
ous NMR/NQR properties. Whether these two phenom- 
ena are related, however, remains an open question. We 
also suggested § that the discrepancy between the ex- 
perimental temperature dependence of 17 W c g/ 89 W and 
the MMP model prediction, which does not include 
quadrupolar relaxation, could be easily explained by as- 
suming an additional quadrupolar relaxation channel for 
plane oxygen. The data of Fig. || support this idea. 

Further information on the nature of 17 W% is pro- 
vided by the pronounced temperature dependence of the 
anisotropy of the effective relaxation rate 17 W e g which 
is not expected within the SSF model Martindale 
et al. H tried to explain this dependence by assuming 
two spin degrees of freedom. In contrast, the uniform 
spin susceptibility monitored by the different plane nu- 
clei scales perfectly thus strongly supporting the view of a 
single spin degree of freedom. Therefore, we believe that, 
below about 200 K, the temperature dependence of the 
17 Wofr anisotropy is caused by the additional quadrupo- 
lar relaxation channel. If so, this enables us to determine 
the anisotropy of W2 ■ 

The magnitude of the W2 rate itself yields an im- 
portant information on the nature of the quasiparti- 
cles. Below 200 K, W2 is much too large to originate 
from simple quadrupolar relaxation due to electron-like 
quasiparticles [^0|. Their fluctuation spectral density is 
smeared out up to frequencies that are extremely high 
as compared to the nuclear Larmor frequency; there- 
fore, the charge fluctuation amplitude at Larmor fre- 
quency, and consequently W2, is very small. The sub- 
stantial quadrupolar relaxation we observe indicates that 
it has to arise either from strongly correlated quasiparti- 
cles as proposed by different theoretical models [^l], or 
from very strong electron-phonon interactions that could 
lead to heavy polaronic-like quasiparticles (^2|. In both 
cases, the slown-down quasiparticle dynamics could pro- 
duce enough low-frequency spectral density necessary to 
account for the observed quadrupole relaxation. 

One nevertheless would expect this relaxation to be 
equaly operative at oxygen and copper sites and, due 
to the larger copper nuclear quadrupole moment, even 
more effective at the later sites. This is not what one 
observes. Whereas there is a substantial quadrupolar 
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contribution to the oxygen relaxation below 200 K, as 
shown in this work, no such contribution has been de- 
tected in the copper relaxation. Since the squared ratio 
of the quadrupole moments of 63 Cu and 17 O is about 60, 
one expects that such a contribution would be easily de- 
tected. We have measured very accurately (by nuclear 
quadrupole resonance) the ratio R of the 65 Cu and 63 Cu 
relaxation rates at the plane copper site in our sample 
at 100 K. Our result, R = 1.1497(35), is in the small 
error limits equal to the squared ratio of the copper gy- 
romagnetic ratios. The high precision of R sets the limit 
for any possible quadrupolar contribution to the Cu re- 
laxation to less than 1 % which is far from what one is 
expecting. 

Finally, let us return to the remark, made at the be- 
ginning of this Letter, that antiferromagnetic fluctuations 
are very well filtered out at the plane oxygen sites. This is 
true if the fluctuations are commensurate with the crys- 
tal structure. However, recent inelastic neutron scatter- 
ing experiments revealed that there are incommensurate 
antiferromagnetic correlations developing in underdoped 
YBa2Cu307_,5 pjj. Most likely such correlations are 
present also in YBa2Cu40s, so one expects much less 
efficient filtering of fluctuations at the oxygen site and a 
more copper-like oxygen magnetic relaxation in contrast 
to what is observed experimentally The discrepancy in 
YBa 2 Cu40s is even more pronounced, since, so far, the 
effective rate 17 W c g was used to analyze the data and 
thus the magnetic relaxation rate W got overvalued. The 
reconciliation of these NMR findings with those gained by 
neutron scattering experiments poses another challenge 
to the theory. 

In conclusion, we have shown that the spin-lattice re- 
laxation of plane oxygen in YBa2Cu40s below approxi- 
mately 200 K is not driven only by magnetic but also by 
quadrupolar fluctuations (i.e. low- frequency charge fluc- 
tuations). In the superconducting state, this newly es- 
tablished quadrupolar relaxation diminishes faster than 
the magnetic one indicating that the underlying relax- 
ation process is strongly influenced by the superconduct- 
ing transition. There are two degrees of freedom involved 
in the low-energy excitations of the electronic system, 
one of them is the single-spin degree, implying that the 
single-spin fluid model is partially correct, whereas the 
other one is the charge degree of freedom with predom- 
inantly oxygen character, since it is not observed at the 
copper sites. 

The partial support of this work by the Swiss National 
Science Foundation is gratefully acknowledged. 
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